%Beta Robust Weight estimation
function [beta]  = BetaRobustWeight(rStock,rIndex,weightfunction)
    beta = robustfit(rIndex,rStock,weightfunction);
    beta = beta(2,:);
end